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Abstract 

Using data collected by the HERA-_B experiment, we have measured the fraction of J/if’s produced 
via radiative Xc decays in interactions of 920 GeV protons with carbon and titanium targets. We 
obtained = 0.32±0.06stat±0.04sys for the fraction of .J/ip from Xc decays averaged over proton- 
carbon and proton-titanium collisions. This result is in agreement with previous measurements and 
is compared with theoretical predictions. 

PACS: [13.85.Ni], [13.85.Qk], [24.85.-bp], [24.85.Eqp] 
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1 Introduction 


The mechanism by which quarkonium states are produced in hadronic collisions is not understood and 
is a subject of current interest. At present, several models exist. The Colour Singlet Model (CSM) 
[H requires that the qq pair be produced in a colour singlet state with the quantum numbers of the 
final meson. The Non-Relativistic QCD factorisation approach (NRQCD) [Q, |3| assumes that a colour 
singlet or colour octet quark pair evolves towards the final bound state via exchange of soft gluons. The 
nonperturbative part of the process is described by long distance matrix elements which are extracted 
from data. Finally, the Colour Evaporation Model (CEM) 1^, assumes the exchange of many soft gluons 
during the formation process such that the final meson carries no information about the production 
process of the qq pair. 

Charmonium production is an attractive test case as the quarks are heavy enough for perturbative 
calculations of the qq production process, yet the cross sections are large enough to be measured with 
good statistics. The dependence of the ratio of production cross sections for different states, e.g. the 
ratio of XcQ and J/'0 production cross sections a{xc)/o'{J/'4’)^ oa ^/s or the projectile allows one to 
distinguish among different models. From the experimental point of view, the specific decay Xc ^ 
is advantageous since the decay signature J/?/) —> (i = can be used as trigger requirement. 

Furthermore, several systematic errors cancel in the ratio, and the only significant difference in the 
detection of the Xc and the J/'0 is the photon reconstruction. Due to the small branching ratio of 
Xco J/^1, (6.6 ± 1.8) • 10-3 @, the XcO contribution to the reconstructed Xc signal can be neglected. 
The Xci and Xc 2 states, with radiative branching ratios of 0.273±0.016 and 0.135±0.011 Q, respectively, 
are separated by 46 MeV/c^. In most experiments the energy resolution is insufficient to resolve these 
two states, so that one usually quotes the ratio 


^xc — 


E cr{Xci)Br{xc 

i=l 






( 1 ) 


Here, is the sum of production cross sections for direct J/V^’s and J/ijj's produced in decays 

of Xc and ip'. In the same way, cr{xci) includes direct Xc production and the feed-down from the p'. 
Contributions from and heavier charmonia are neglected. 

While this “inclusive” ratio is usually quoted in the literature, one can define the ratio for direct 
Xc production over direct J/p production 


E '^{Xci)dirB'f(^Xci 

jydir _ 2=1 _ 

(^{JlP)dir 


J/Px) 


( 2 ) 


can be derived from and the known ratio of p' to J/p production cross sections ||^] and known 
branching ratios [^. 

The experimental situation is unclear, and the uncertainties are large particularly for proton induced 
reactions where the few existing measurements of |^] differ strongly. Measurements made with 
pion beams have higher precision and may indicate that R^^ increases with ^/s. For photon and 
electron-induced reactions, only upper limits for have been reported |1C|. 

We report here a new determination of Ry-^ in interactions of 920 GeV protons with carbon and 
titanium nuclei. The Xc is observed in the decay Xc ^ J/Pl using the value AM, 

which is the difference between the invariant mass of the system and the invariant mass of the 


^In the following, the notation “Xc” indicates the sum of the three states XcO, Xci and Xc2- 
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lepton pair £+£ 


AM = Mii+i-j) - M{i+t 


( 3 ) 


Here, the uncertainty in the determination of the J/if} mass essentially cancels. An excess of events with 
respect to the combinatorial background determines the number of Xc candidates from which the 


inclusive” value, can be calculated as follows: 


Rx. - ^ 


Ny 


J/ij) ‘ ^7 


* Pei 


( 4 ) 


where N is the total number of reconstructed J/ij) ^ decays used for the Xc search. The factor 
is the photon detection efficiency. The value represents the ratio of trigger and reconstruction 


efficiencies for J/ip’’s from Xc decays and for all J/ip's: 


Ps = 


J/P J/p 

trig <^reco 

i-Xc i-Xc 

^trig ^reco 


( 5 ) 


Since the kinematics, triggering and reconstruction of direct J/tp's and J/tp^s from Xc decays are very 
similar, p^ is close to unity. 


2 Detector, Trigger and Data Sample 

HERA-H is a fixed target experiment operating at the HERA storage ring at DESY. Charmonium and 
other heavy flavour states are produced in inelastic collisions by inserting wire targets into the halo of 
the 920 GeV proton beam circulating in HERA. The pN (N = p,n) c.m.s. energy is ^/s = 41.6 GeV. 

The detector is a magnetic spectrometer emphasising vertexing, tracking and particle identification, 
with a dedicated J/ip-tiiggeT. The components of the HERA-R detector used for this analysis include 
a silicon strip vertex detector (VDS), honeycomb drift chambers (OTR), a large acceptance 2.13 T-m 
magnet, a finely segmented “shashlik” electromagnetic calorimeter (EGAL), and a muon system (MUON) 
consisting of wire chambers interleaved with iron shielding which detects muons with momenta larger 
than 5 GeV/c. The EGAL is divided into three radial parts with decreasing granularities, of which two, 
the “inner” and “middle” sections, are used for the measurement described here. The performance of 
the calorimeter is described in Ref. |T^]. The HERA-R detector allows an efficient reconstruction of 
particles with momenta larger than 1 GeV/c, including y’s and tt^^’s, within the acceptance. A detailed 


description of the apparatus is given in Ref. 

The HERA-R target station houses 8 target wires which can be moved independently into the beam 
halo. Their positions are steered such that the proton interaction rates are equalised for the targets 
in use. The data presented here were obtained using a carbon wire and a titanium wire separated by 
3.3 cm along the beam direction. The resolution of the reconstructed dilepton vertices of 0.6 mm along 


the beam direction |13] allows a clear association of the interaction to a specific target wire. The analysis 
presented is based on data collected during a short commissioning run in summer 2000. About half of 
the data was taken with a single carbon wire; the second half was taken with carbon and titanium wires 
together. The proton-nucleus interaction rate was approximately 5 MHz. 

The trigger selects and e^e“ pairs, the latter with an invariant mass larger than 2 GeV/c^. For 

a muon candidate the trigger requires at least 3 MUON hits in coincidence with 9 OTR hits consistent 
with a particle track with a transverse momentum between 0.7 GeV/c and 2.5 GeV/c. The electron 
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trigger requires that the transverse energy deposited in the calorimeter^ by the electron candidates be 
greater than 1.0 GeV and that at least 9 OTR hits confirm the track hypothesis. Both muon and electron 
candidates have to be confirmed by a track segment in the vertex detector with at least 6 hits. For the 
data described here, the trigger acceptance for J/f/^’s was limited to the xp range —0.25 < xp < 0.15, 
xp being the Feynman’s x variable. For more details concerning the trigger and the data sample of the 
year 2000 run, see Ref. Q. 

The data is divided into four separate samples: or e+e“ final states, each originating from 

either carbon or titanium target wires {fi-C, e-C, /r-Ti and e-Ti). 


3 Monte Carlo Simulation 


At present, NRQCD is the favoured approach to describe charmonium formation. It is therefore used to 
generate our signal Monte Carlo sample. To estimate the model dependence systematics we use the CSM. 
Since the GEM does not make any conclusive predictions for the differential charmonium production 
cross sections, we have not used it in the simulations. 

The Monte Carlo simulation (MG) of events is done in three steps. First, a J/'0 or Xc is generated 
using PYTHIA 5.7 In the simulations, the CTEQ2L parton density function [jl^ ] and the c quark 
mass rric = 1.48 GeV/c^ are used. The sum of the transverse momenta of the reaction products must 
exceed 0.5 GeV/c. Any polarisation is neglected. Eor NRQCD, the differential cross sections and long 
distance matrix elements are taken from Ref. [^. Eor CSM, the differential cross sections are taken 
from Ref |l|. During the second step, the energy remaining after the charmonium generation is used to 


simulate the rest of the pA interaction using FRITIOE 7.02 [16|. Finally, the J/ip event is combined 
with n other inelastic interactions to simulate several interactions per event, as observed in the data. 
The number n is distributed according to Poisson statistics with a mean value of 0.5 determined from 
the mean experimental interaction rate. 


The detector response is simulated using GEANT 3.21 l]!^ and includes the measured hit resolution, 
the mapping of inefficient channels, and electronic noise. The simulated events are processed by the 
same trigger and reconstruction codes as the data. The simulation has been checked to ensure that it 
accurately describes the detector, both in terms of geometric acceptance and material composition (see 
sect. |4.6| ). From the MG we expect a mass resolution for the Xc signal of 45 MeV/c^, which is insufficient 
to separate the Xci and Xc 2 states. 

MG studies show that the trigger and reconstruction efficiencies for J/ijj ^ are indeed similar 
for both direct J/V^’s and for those originating from Xc ^ J/''Pi decays. We obtain = 0.95 ± 0.02 for 
the NRQGD and ps = 0.97±0.01 for the GSM. For the measurement we use the NRQCD value pe = 0.95 
and consider the difference between the two values as a measure of the systematic uncertainty of p^ (see 
sect. 

The Monte Carlo sample used in the analysis is about six times larger than the data sample. 


■^The transverse energy is defined as cluster energy multiplied by the transverse distance to the beam axis and divided 
by the cluster-target distance. 


6 





4 Data Analysis 


4.1 Method and Selection Criteria 

The analysis consists of the reconstruction of the J/ij} events, the search for the photon candidates in 
the ECAL, and the determination of the invariant mass of the J/'0 and photon candidates within the 
event. The selection criteria for the Xc are tuned to maximise the quantity Nji^ . Here, Njj^ is 

the number oi Jjil) candidates above background found in the data, and is the photon reconstruction 
efficiency determined from MC simulations. N is the number of all events, including background, found 
in the measured AM distribution within a window of two standard deviations (determined from MC) 
around the expected Xc position. The procedure is applied for all cuts described below. 

4.2 Selection 

In the offline analysis, a track is selected as a muon candidate if its transverse momentum is greater than 
0.7 GeV/c and the muon likelihood, derived from the ratio of the expected and found MUON hits, is 
greater than 0.001. The latter removes hit combinatorics which satisfy the trigger while keeping nearly 
all good muons. 

A track is identified as an electron candidate if (a) the transverse energy {Et) is greater than 
1 GeV; and (b) it has \E/p — 1| < 0.3, where E is the energy deposited in the calorimeter and p 
is the track momentum. The cut on E/p corresponds to about 3.3 standard deviations of the electron 
E/p distribution. To further reject the background from hadrons, a search for associated bremsstrahhmg 
photons emitted in the region upstream of, or inside, the magnet is performed for each electron candidate. 
Thus, an isolated electromagnetic cluster is required in the area where the bremsstrahlung would hit the 
EGAL. The energy of the bremsstrahlung cluster is added to the energy of the electron candidate. The 
requirement of an associated bremsstrahlung photon candidate for each of the two electron candidates 
of the J/i/ e^e~ decay has an efficiency Shrem of about 20% (about 45% per electron) and suppresses 
the background by a factor of 45. These values are obtained by comparing the J/i/ and background 
rates under this requirement with those for the case that at least one of the two electrons is associated 
with a bremsstrahlung cluster, and they are also confirmed by MC studies (see alse sect. 4.6). 

The assignment of the J/ijj candidates to a target wire is based on the position of the reconstructed 
dilepton vertex. The x^ probability of this vertex is required to be larger than 0.005 to eliminate spurious 
events. 

The invariant mass is calculated for each opposite-charge lepton pair. The resulting mass distributions 
are shown in Eig. || for each of the four samples. The signal observed for J/if) ^ events is Gaussian 

while the J/'i/ ^ e'’'e“ signal has an asymmetric bremsstrahlung tail. In both cases the background 
underneath the signal comes mainly from misidentified hadrons and conversions. The background shape 
is either described by an exponential distribution or by an exponential multiplied by a second 

order polynomial distribution (e^e“). The shape is confirmed by the invariant mass distribution of the 
same sign candidates in the muon case, and by fitting the distribution that results from using all trigger 
candidates (mostly hadrons) in the electron case. Only J/i/ candidates within a two standard deviation 
(2 a) window around the J/i/ mass are considered for the analysis. In the electron case, a is taken from 
the high mass part of the signal which is Gaussian. The numbers of J/i/ candidates obtained from the 
fit and corrected for the mass window are shown in Table || for the four samples. 
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Figure 1: Dilepton invariant mass spectrum for each of the four samples (//-C, /r-Ti, e-C, e-Ti). The 
dashed lines show the estimated background under the signal. See text for the details on the fits 
(solid lines). The selection criteria are described in sects. and |4.3| . 


4.3 Particle Multiplicities 

Detector occupancies have a considerable impact on the Xc reconstruction: large calorimeter occupancies 
lead to more combinatorial background. However, detector occupancies are correlated with particle 
multiplicities and thus depend on the underlying physics of the event. Since the cross section for J/'ij) and 
Xc production are of the same order of magnitude, and the kinematic dependence of direct J/ip and those 
from Xc decays are similar, we assume both types of events to have similar particle multiplicities excluding 
the decay products of charmonium. The charged particle multiplicity is proportional to the number of 
tracks reconstructed in the VDS. We eliminate especially “busy” events which tend to contribute more 
to the background than to the signal. Based on the multiplicity distributions we require that there 
are not more than 30 (34) VDS tracks in events with a J/'0 candidate originating from the carbon 
(titanium) target. We also require less than 19 clusters in the ECAL. The upper limit on the number of 
tracks reduces the background under the J/'0 signal while the cut on the number of clusters limits the 
combinatorial background under the Xc- The numbers of passing the multiplicity cuts are given 

in Table |l|. 

4.4 Photon Selection 

Each cluster in the ECAL with Et >0.1 GeV that is not associated with the leptons from the J/'ijj, is 
considered as a photon candidate. The area of the ECAL closest to the proton beam pipe x^/4 + < 

484 cm^ (or equivalently : 0^/4 + 0^ < 265 mrad^) is excluded, as the occupancy in this region is high 
(up to 30%). Hadronic background is reduced by requiring that the ratio of the central cell energy to 
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the total cluster energy (Ecentr/E) be greater than 0.6. In order to suppress background due to soft 
secondary particles and noise clusters, an energy cut E > 3.0 GeV is applied. A charged track veto 
is not applied, due to a 44% probability of the photon to convert in the detector material downstream 
of the magnet. The material of the detector in front of the ECAL causes photon conversion, and thus 
losses of photons from Xc’s. We determine these losses using MC simulations. 

Since the relative momenta of the J/ij} and the photon from the Xc decay are correlated, a cut in 
the acceptance of the J/if: affects the acceptance for the photon as well. The electron sample, with a 
slightly larger acceptance close to the beam as compared to the muon sample, also contains more energetic 
photons than the muon sample. The different samples have different kinematics and acceptances, leading 
to different photon detection efficiencies which are determined for each sample separately using MC 
simulations (see Table ||). The uncertainty in the photon detection efficiency arises mainly from the 
finite MC statistics; however, this uncertainty is insignificant compared to the statistical error on 

4.5 Xc Reconstruction 

The AM distributions for all combinations ol J/il) and photon candidates for the carbon samples are 
shown in Fig. ||. The distributions show a signal corresponding to the sum of the two charmonium states 
Xci and Xc2- 

Possible sources of background are random combinations oi J/tp and photon candidates, decays of 
heavier mesons into J/pX^ and the radiative decay J/'0 ^ e'’'e“ 7 . The fraction of J/'p's originating 
from p' decays is about 8% [^. The fraction of photons arising from p' J/pTr^n^ decays which 
pass the energy cut E > 3.0 GeV is negligible, as is the fraction of pions misidentified as photons 
from p' — > J/p'K~n~^ decays. The fraction of J/p's resulting from decays of S-mesons, T, Xb^ v'ci 
XcO is negligible as well. The photon from the radiative decay J/p ^ e'^e~'y, mostly oriented along 
the direction of one of the leptons, is indistinguishable from bremsstrahlung and, as such, is taken into 
account. In the muon sample, bremsstrahlung clusters are neither expected nor found above background. 
The fraction of such radiative decays is considered to be negligible. Thus the background consists mainly 
of random combinations of J/p and photon candidates. 

The shape of the dominantly combinatorial background in the AM distribution is obtained by com¬ 
bining J/p candidates with photon candidates from different events with similar multiplicity and ap¬ 
plying the standard selection cuts. These “mixed events” reproduce the shape of the AM distribution 
everywhere except in the Xc signal region (see solid line in Fig. left panel). Similar results are obtained 
when events in the sidebands of the J/p mass region are combined with photon candidates. Since the 
experimental resolution is of the same order as the mass difference between Xci and Xc 2 states and the 
statistics is limited, we use a single Gaussian to describe the signal. In the fit, the position and normal¬ 
isation of the Gaussian, as well as the normalisation of the background, are left free. The width of the 
Gaussian is fixed to the value predicted by MC based on the NRQCD approach (45 MeV/c^), where the 
production cross section ratio of Xci and Xc 2 is approximately 0.65. The position of the Gaussian agrees 
well with the value expected from MC. The background normalisation is also treated as a free parameter 
when we fit the number of Xc J/Pl decays. The background subtracted distributions are shown in 
the right panel of Fig. |2|. The significance of the signals seen in the /x-C and e-C samples is about three 
standard deviations. The obtained number of Xc events as well as the number oi J/p events and the 
photon detection efficiency are summarised in Table ||. 

Taking into account the high background level and the ratio of to X observed in the carbon 
samples, we do not expect to see a signihcant Xc signal in the smaller titanium samples. The results 
obtained for the titanium sample with the same procedure are shown in Fig. ^ The starting value for 


9 



Figure 2: The AM = M{i^£~'y) — M{£~^i~) distribution for samples /r-C and e-C. In the left-most plots, 
the points represent data and the solid lines represent the combinatorial background estimated by event 
mixing. The right-most plots show the signal after background subtraction. See text and Table for 
the details on the fit. 


the peak position has been taken from the fit of the carbon data. Although the signals are marginal, the 
values obtained from them are compatible with estimates from the carbon samples (see Table |^. 
As a cross check, all four samples are combined together as shown in Fig. ^ The value of = 
380 ± 74 obtained from this distribution agrees within the errors with the sum of the values obtained 
from the individual samples. 

4.6 Study of the Systematic Uncertainties 

The systematic uncertainty in the yield of J/V' —> e'’'e“ candidates due to the background description 
is estimated to be 5%, whereas the uncertainty is negligible in the muon case. 

The model dependence of the relative efficiency pe for all J/tp’s and J/'ip's from Xc has been studied. 
A 2% difference of is found for the two models (NRQCD and CSM) for each of the J/ip leptonic 
decay modes. For the same models a difference of the photon detection efficiency integrated over 
all Xc states is found to be 4%. In both cases the observed difference is treated as an estimate of the 
corresponding systematic uncertainty. The overall systematic error accounting for the model dependence 
of the selection efficiency is 5%. 

To confirm the MC description of the detector material composition and acceptance which affects 
the photon detection efficiency, we compare the bremsstrahlung tag probability y/etrem determined from 
the data and MC. The values obtained, 0.44 ± 0.02 and 0.43 ± 0.01 for the data and MC, respectively, 
are compatible within one standard deviation. The systematic uncertainty due to photon losses is 
conservatively taken to be 2%. 
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Figure 3: Same as in Fig. for the ^-Ti and e-Ti samples. 
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Figure 4: Same as in Fig. 0 for all data combined. 
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Table 1: The number of selected J/ip events the number of J/ip's passing the multiplicity cut 

{Nj/^), the number of Xc’s observed pe^' degree of freedom for the AM ht, photon detection 

efficiency (e-y), and the result for for each of the four event samples. The quoted error on is 
statistical, excluding the systematic uncertainty in e^. 



fi-C 

e-C 

/i-Ti 

e-Ti 

j\jselected 

1760 ± 48 

1380 ± 69 

765 ± 31 

512 ±41 

^j/i) 

1510 ±44 

1180 ± 59 

643 ± 29 

382 ± 32 

Nx. 

159 ± 47 

121 ± 38 

59 ±33 

31 ± 27 

X^ln.d.f 

28/35 

34/35 

27/35 

48/35 

Ey (%) 

27.3 ± 1.1 

32.8 ± 1.5 

24.4 ± 1.8 

32.7 ±2.6 

^Xc 

0.37 ±0.11 

0.30 ± 0.09 

0.36 ± 0.20 

0.23 ±0.21 


The effect of systematic uncertainties of the ECAL calibration is studied by using MC simulations. 
The level of possible uncertainties is determined from the data using the vr^ signal. The uncertainty of 
the TT^ calibration is used then in the MC simulations to determine its effect on the detection of the Xc- 
The systematic uncertainty on Ry^ due to this effect is estimated to be 1%. 

Correlated electronic noise in the calorimeter can shift and widen the AM distribution of the Xc 
signal. A cluster reconstruction procedure based on the known correlation between channels is developed 
to compensate for this effect. The numbers of events observed with and without this algorithm in the 
data agree with each other within the statistical errors. MC studies show that the corresponding relative 
systematic uncertainty is 3%. 

The width of the Xc signal in the AM distribution depends on the ratio of Xci and Xc 2 and on the 
detector resolution, mainly that of the ECAL. When the width is left free to vary in the carbon data fit, 
the resulting width agrees within one standard deviation with the nominal one. A systematic error on 
Ry^ of 6% is assigned based on an MC study of the signal resolution dependence on the Xci to Xc 2 ratio 
and the ECAL resolution. 

The stability of our results with respect to variations in the selection criteria is studied separately 
for the different samples. The ratio Ry^ is measured as a function of the cuts on VDS track multiplicity, 
ECAL cluster multiplicity, photon energy E, and ratio Ecentr/E. The variation of the cut on the photon 
energy E results in a variation of Ry^ of 6%, which is taken as an estimate for the systematic uncertainty. 
The dependence on other cuts is negligible. 

The systematic uncertainty on due to the finite MC statistics is 3% as estimated from the weighted 
average of the values in Table |I[ 

The polarisation of the Xc affects the reconstruction efficiency of Xo however with the present statis¬ 
tics we are not able to determine it. Thus, like the previous experiments, we have assumed no polarisation 
and have neglected the uncertainty related to it. 

Assuming that all individual systematic errors are uncorrelated, an estimate of the total systematic 
uncertainty on Ry^ is 11% (see Table ^). 
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Table 2; Contributions to the relative systematic uncertainty. For the calculation of the total uncertainty 
the correlation in the systematic errors of the different samples is taken into account. 


Uncertainty 

J/i!) background shape (e-C, e-Ti only) 

5% 

Model dependence 

5% 

Photon losses 

2% 

ECAL calibration 

1% 

Correlated noise in ECAL 

3% 

AM resolution 

6% 

Dependence on cuts 

6% 

Finite MC statistics {e^) 

3% 

Total 

11% 


5 Results 

The values of Ry.^ obtained for all four samples are listed in Table The results for the two carbon 
samples agree with each other within the statistical errors. The results obtained from the titanium data 
are consistent with those obtained from the carbon data. 

Although nuclear dependence effects might be present in at the few percent level for the targets 
used here ||T^ , they are beyond the statistical accuracy of the present measurement. We therefore average 
the results for the four samples obtaining: 

{R^J = 0.32 ±0.06,tet± 0.04,2,,. (6) 

The hrst uncertainty listed is statistical only, whereas the second uncertainty is systematic. 

In order to extract the ratio of the “direct” Xc and J/il: production, we use 

ndir _ _ 1 — R^Bl _^ 

1 - + R^B2 ’ ^ ’ 

where R^ = a{'ijj')/a{J/'tp) = 0.094 ± 0.035 is taken from Ref. 0 and corrected for the branching ratios 

2 

[^. Bi is the branching ratio Br{'ip' —> J/tpX) and B 2 is the sum of branching ratios Br{'ip' 

Xca)Br{xci [i- We obtain R^l = 0.50 ± 0.15,tai ± 0.08,^,. 

Our result for Ry^^ (eq. (6)) is compatible with most of the previous data (|^, ^ and |^), as shown 
in Table |3| and Fig. 0 Due to the relatively large uncertainties, especially for the data on proton induced 
reactions, a flat energy dependence, as predicted by GEM 0, cannot be ruled out. Similarly, the slope 
of the energy dependence as predicted by the Monte Carlo based on NRQCD (see sect. 0 is also 
compatible with the data. However, the predictions of NRQCD seem to fall below the data, which might 
indicate that the NRQCD long distance matrix elements for Xc production Q used for the calculations 
are underestimated. On the other hand, CSM predicts values for R^^ which are larger than most of 
the data. More precise measurements, especially of proton induced reactions are needed to conclusively 
discriminate among these models. 
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Figure 5: Comparison of our measurement of (closed triangle) with those of other pp, pA Q 
(closed circles) and vrp, ttA (open circles) experiments. The CDF result ||^ is not shown, since its 
kinematic acceptance differs strongly from the other experiments. The value quoted for exp. E771 has 
been calculated from the published cross sections and branching ratios [Q. The error bars include 
statistical and systematic uncertainties. Also shown are predictions for pN and ttN interactions obtained 
from Monte Carlo based on the NRQCD (solid), CSM (dashed) (see sect. The CEM Q, |5| 
predicts a constant value. The dot-dashed line is the average of all measurements. 
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Table 3: Previous -kA |^], pA P, pp |^| and HERA-S measurements of the value. The value quoted 
for exp. E771 has been calculated from the published cross sections Q and branching ratios 


Exp. 

coll. 

(GeV) 

^Xc 

IHEP140 

7r~p 

8.5 

0.44 ±0.16 

WAll 

7r“Be 

18.7 

0.30 ±0.05 

E610 

7r“Be 

18.9 

0.31 ±0.10 

E673 

7r”H2, 7r“Be 

20.2 

0.70 ±0.28 

E369 

7r“Be 

20.6 

0.37 ±0.09 

E705 

7r“Li 

23.8 

0.37 ±0.03 

E705 

vr^Li 

23.8 

0.40 ±0.04 

E672/706 

7r“Be 

31.1 

0.443 ± 0.041 ± 0.035 

E610 

pBe 

19.4, 21.7 

0.47 ±0.23 

E705 

phi 

23.8 

0.30 ±0.04 

E771 

pSi 

38.8 

0.74 ±0.17 

R702 

pp 

52, 63 


ISR 

pp 

62 

0.47 ±0.08 

CDF 

pp 

1800 

0.297 ±0.017 ±0.057 

HERA-R 

pC, pTi 

41.6 

0.32 ± 0.06 ± 0.04 


6 Conclusions 

A measurement of the ratio of J/'i/) produced via radiative Xc decays to all produced J/V’ allows one 
to quantitatively test different models of quarkonium production. We present a new result from the 
HERA-R experiment for the fraction of J/ip's originating from radiative decays of Xci Xc 2 states 
produced in pC and pAi interactions. The fraction of J/ip's in the range of —0.25 < xp < 0.15 originating 
from radiative Xc decays is determined to be = 0.32 ± O.OG^tat ± O.Od^j^s, and consequently, the ratio 
of the cross section of directly produced Xc’s decaying into J/ip to the cross section of directly produced 
J/'0’s is = 0.50 ± 0.15stat ± O.OSsys in the above xp range. The result has been obtained with C 
and Ti targets and detecting the J/'0 decay modes into electrons and muons. Our result for R^^ agrees 
with most previous proton and pion beam measurements, neglecting any possible energy dependence. R 
agrees also with the predictions of the Non-Relativistic QCD factorisation approach (NRQCD), whereas 
it falls significantly below the predictions of the Color Singlet Model (CSM). 
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